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THE INERTIA COEFFICIENTS OF AN AIRSHIP IN A FRICTIONLESS FLUID. 

By H. Bateman. 


SUMMARY. 

The following investigation of the apparent inertia of an airship hull was made at the 
request of the National Advisory Committee for Aeronautics. The exact solution of the aero- 
dynamical problem has been studied for hulls of various shapes and special attention has been 
given to the case of an ellipsoidal hull. In order that the results for this last case may be 
readily adapted to other cases, they are expressed in terms of the area and perimeter of the 
largest cross section perpendicular to the direction of motion by means of a formula involving 
a coefficient K which varies only slowly when the shape of the hull is changed, being 0MS7 
for a circular or elliptic disk, 0,5 for a sphere, and about 0.25 for a spheroid of fineness ratio 7. 
For rough purposes it is sufficient to employ the coefficients, originally found for ellipsoids, for 
hulls otherwise shaped. Wheit more exact values of the inertia are needed, estimates may be 
based on a study of the way in which K varies with different characteristics and for such a 
study the new coefficient possesses some advantages ov'er one which is defined with reference to 
the volume of fluid displaced. 

The case of rotation of an airship hull has been investigated also and a coefficient has 
been defined with the same advantages as the corresponding coefficient for rectilinear motion. 

I. INTRODUCTION. 

It follows from Green’s analysis that when an ellipsoidal body moves in an infinite incom- 
pressible inviscid fluid in such a way that the flow is everywhere of the irrotational, continuous 
Eulerian type, the kinetic energy of the fluid produces an apparent increase in the mass and 
moments of inertia of the body. The terms mass and moment of inertia are used here in a 
generalized sense because it appears that the apparent mass is generally different for different 
directions of motion and the apparent moment of inertia different for different axes of spins 
For this reason it seems better to speak of inertia coefficients, these being the constant coeffi- 
cients in the expression for the kinetic energy in terms of the component linear and angular 
velocities relative to axes fixed in the body. 

The idea of inertia coefficients may be extended to bodies of any shape and to cases in 
which there is more than one body or in which the fluid is limited by a boundary. Generalized 
coefficients may be defined, too, for cases in which there is circulation round some of the bodies 
or boundaries and values can eventually be obtained which should correspond closely to the 
values of the inertia coefficients for the motion of a body in a viscous fluid. 

The inertia coeflTicients of airship hulls are useful for the interpretation of running tests 
and in fact for a dynamical study of any type of motion of an airship, whether steady or unsteady. 
The coefficients are needed, for instance, in the study of the stability of an airship by the method 
of small oscillations * and for a computation of the resulting momenta in various types of steady 
motion. 

For the case of motion of translation with velocity V the kinetic energy, T, of the fluid 
is usually expressed in the form 

T=\lcmV^ 

where m is the mass of the fluid displaced by the body and ^ is a numerical coefficient whose 
value is known in certain cases. A value of k for an airship hull is generally found by choos- 

* For the literature on this subject reference may be made to a paper by R. Jones and D. H. Williams, British Aeronautical Research Com- 
mittee, R. M. 751. June, 1921. 
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ing an ellipsoid with nearly the same form as the hull and calculating the value of k for the 
ellipsoid. This method is to some extent unsatisfactory because the coefficient k varies con- 
siderably with the shape, being infinite for a circular disk, 0,5 for a sphere, and 0,045 for a 
prolate spheroid of fineness ratio 6. For this reason an alternative method is proposed in which 
the kinetic energy of the fluid is expressed in terms of quantities relating to the master section 
of the hull by means of a formula involving a numerical coefficient K which varies only slowly 
with other characteristics such as the fineness ratio. The proposed expression is 




where S denotes the area and I the perimeter of the greatest cross section of the hull by a plane 
perpendicular to the direction of motion; p is the density of the fluid and K the new coefficient 
which is apparently greatest for a circular or elliptic disk. 

In the case of a spheroid moving in the direction of its axis of symmetry the way in which 
k and Z vary with the fineness ratio is shown in Figure 1. In Figure 2 the corresponding curves 
have been drawn for a hull bounded by portions of two spheres cutting each other orthogonally. 
The high value of K when the two spheres are equal is undoubtedly caused by the presence of 



the narrow waist, while the sudden drop in value indicates the effect of a lack of fore and aft 
symmetry. The curves for Zhave an advantage over those for k in indicating more clearly the 
effect of a change in shape. The effect of a flattening of the nose of the hull has been studied 
by considering the case of a surface of revolution whose meridian curve is a limacon. The effect 
is only slight, as is seen from the table in Section IV. In the case of an airship hull spinning about 
a central axis in a plane of symmetry the kinetic energy can also be expressed in terms of general 
characteristics by using a formula involving a coefficient /iT, which varies only slowly with 
the shape. The proposed formula is 



K',R^ 


is,-s,y 

~~ i 


wdiere is the angular velocity about the axis of spin, which we take as the axis of x, is the 
maximum radius of gyration of a meridian section about the axis of x, Sy and Sz are the areas of 
central sections perpendicular to the axes of y and z and I is the perimeter of the meridian section 
with the greatest perimeter, a meridian section being cut out by a plane through the axis of spin. 

This formula has been constructed from the known formula for an ellipsoid with the axes 
of coordinates as principal axes. To adapt it to a hull of a different shape a suitable set of axes 
must be chosen. The principal axes of inertia at the center of gravity may, perhaps, be used 
with advantage. 
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The coefficients K and K' will now be computed in some cases in whiUh the aerodynamical 
problem is soluble. In particular they will be computed for the following cases: 

(1) Disk moving axially. 

(2) Prolate spheroid moving longitudinally. 

(3) Prolate spheroid moving laterally. 

(4) Oblate spheroid moving in the direction of its axis of symmetry. 

(5) Oblate spheroid moving at right angles to its axis of s^mimetry. 

(6) Solid formed by two orthogonal spheres. 

(7) Solid formed by the revolution of a lima^on about is axis of symmetry. 


II, THE INERTIA COEFFICIENTS FOR AN ELLIPSOID. 


When the viscosity of the fluid is neglected and the motion is treated as irrotational there 
is no scale effect. This means that if we increase the velocity of the body in the ratio s:l , keeping 
its size constant, the velocity at any point of the fluid changes in the same proportion. A 
similar remark applies to the case in which the body is spinning about an axis instead of moving 
with a simple motion of translation and in the more general case in which a body lias motions 
of both translation and rotation the kinetic energy, T, can be expressed in the form^ 


2 T= Au^ + Bv^ + CiD^ + 2A'vw + 2B'wu + 2 C^uv + Pp^ + -1 Rr^ + 2P'qr + 2Q'rp + 2R'pq 

+ 2p {Fu^-Gv-\-Hwy^-2q {F'u-{-G'vA HGd) +2r {F"u-\-G''v + R''v)) , 

where ('U, v, w) are the component velocities of a point fixed in the body and (p, q^ r) are the 
angular velocities of the body about axes through this point that are likewise fixed in the body. 
The coefficients A, /?, C, A', B\ O', P, Q, R, P', Q', R', F, G, H, F', F", G' , IP' are 

constants which are called the inertia-coe ficlents of the body relative to these axes. This expres- 
sion for the kinetic energy has been used also in cases in which the velocities are variable and the 
determination of the inertia coefficients is evidently a matter of some importance. 

The inertia coefficients are usually found by writing down the velocity potential or stream- 
line function which specifies the flow and calculating the kinetic energy by means of an integral 
of type 


over the surface of the body, being the velocity potential, p tlic density of the fluid ami dn 
denotes an element of the normal to the surface dS drawn into the fluid. A different integral 
may be used when the stream-line function is known, but in many cases integration is unneces- 
sary, for Munk^ has remarked that in the case of a simple velocity of translation the fluid motion 
may be supposed to arise from a series of doublets and that the sum of the moments of all these 
doublets has a component in tlie direction of motion which is proportional to the sum of the 
kinetic energy of the fluid and the kinetic energy which the fluid displaced would have if it moved 
like a rigid body with the same velocity as the body. The sum of the masses of the fluid and 
the fluid displaced has been called the complete mass. 

The inertia-coefficients are well known for the case of an ellipsoid with semiaxes a, h, c wlien 
the axes of reference are the principal axes of the ellipsoid. We have in fact^ 


“0 etc 


m=l^pahc,A' = B'^C'^P’ = Q'^R' = F^F'=F"=G = G' = G"^H=n'^II" = o. 

O 


where 


X oo ^00 ^ 

= (6’ 4) A’ = 


d\ 

(c* + X)A 


, A = [(ft» + X)(6’ + X)(c* + X)]!. 


J See Lamb’s hydrodynamics. 

» Technical N^e No. 104, National Advisory Committee for Aeronautics. July, 1»22. 
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The complete cbefficients of inertia for motion of translation are 

A* = ^— , «* = 

2-«o 

The coefficient Ic defined by the equation 


.4* = 2— . = ^ = S-- 

2 — Uq 2 — Pq ^ To 


7.- 0^0 

2-a, 

has been tabulated by Professor Lamb in a number of cases.’ We have extended his tables 
and have also tabulated the coefficient K defined in I. The different special cases of an ellipsoid 
will now be discussed. 

1. Elliptic disk . — In this case k is infinite but the kinetic energy is finite. To find an 
expression for K we write 

oo 

dX _ X ^ 1 ■ 3 • • • {2n — xj _ , j, 

(a*+X)A“^ 2A-?- 

n=o 

2 




X 


abc dX 

(a^+xp(c* + X)" 


{■ 


d\ 


^ - tan~‘ - 7 ^- - 

-aL2 c 


2a TT 


(a^ +X)*(c^ +X) 

when a is small. Differentiating once with respect to a and n times with respect to we get 


(-l)“n! 


Hence 


I 


ad\ 


+X)’ (c^ +X) 


n+l /.2n+2 


V ^2n+3 \ ^ 


2“ 


2br.^l , L3 /C--6JY 
•~T L "^2 c * "^ 2 - 4 ^ / 




1 irabr l-3c^-6* l-y-o/d-h^V 

■J c* ^ 2^ 2^-4^ \ cW 

(c^cos^g + &-sin^>) i of a] =2 [l - 


+ 


dd 


] 


approximately, where I is the perimeter of the ellipse with semiaxes b and c. Hence finally we 
obtain 


and* 


2T=~^^, S = icbc 

OTT I 


^=^1^0.637 

TT 


The distribution of doublets may be found from the well-known expression for the potential 
We have for an ellipsoid 

abc xU I d\ 




As a~^0 we have 


^X 


y" +■ 


4> = — - X U 


2-«o I («^+X)a' a^ + Xo b^ c^-fX( 

dX 


= 1 


^oo 

^ I x*w 

t/x. 


J.2 y2 2^ 

+x)*(?-fX)*’ x/ft^+Xo'^B^+T," * 


Putting \ = x^s and making x^O we find that the value of <f> on the disk is 


« British Advisory Committee for Aeronautics,' R. & M. No. 623, October (1918). 

» In the case ol an inanitely long strip bounded by two parallel lines the value of K is 0.589. 
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This must be equal to the 2ir times the moment per unit area of the doublets in the neighbor- 
hood of the point (0, y, z) of the disk. Hence the expression for the potential is equivalent to 


26c , f f [* 5] 

irl J [x^ + iy-yj^'+iz-zor]- 

and this formula shows the way in which the potential arises from the doublets. The complete 
energy is in this case 


^-3;" [^t^! J[‘- I’-S] 

in accordance with Munk’s theorem. To verify this result we put 

= hs cos 0 ), — cs sin w 

then 

J* 1 - dy^dz^ ^ Vl — dec = he. 

With the above substitution the expression for the potential may be written in the form 

2hc C f - , r^do^ 

^ J ^ j 

iy — hs cos + (2 — cs sin oj)^ 

and may be compared with the corresponding expression for the oblate spheroid. For the case 
of the circular disk (h = c) the stream-line function may be obtained by replacing z in the above 
formula by — (y^-{-z^). When an elliptic disk spins about the axis of y the kinetic energy is 
given by 


2T= ~ TTpe^ti 
1 0 ^ 


r "-J1. 

J CO 


(1 +cos^ B)dd 
cos^ 6 + h^ sin^ 0)^ 


where fiy is the angular velocity. In the case of a circular disk the kinetic energy is 


The coefficient K' thus has the value 


^t = :^ = 0.318. 

7T 


2. Prolate spheroid . — In the case of a prolate spheroid moving in the direction (»f its axis of 
symmetry, we have (Lamb, loc. cit.) 

2(1 \+e 1 

«o= ^3 

where e is the eccentricitv of the meridian section and so 


The velocity potential is 


h = c^a^J\ — < 


where 


_ ah^ Ux r 


d\ 

(a’ + X)* (6* -, X) 


J.2 y^ + z\ 

+ \ V + \ 
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By introducing the spheroidal coordinates 

x = li CO cos Wj 2 = CO sin co = h {1 1)’o h = ae 

we may write this in the form 


^diere 


4 ,, = A P, (m) <^.(f ) = A M [2 f ^ - 1 ] 


The velocity potential may also be expressed as a delinite integral 

<t>A 


H 


sds 

[ {x-hsy + y-+z-\^ 


w'hich indicates the way in which it may be imagined to arise from a row of sources and sinks 
on the line joining the foci. This result may be obtained by writing 

'+1 

f(.s)ds 


Mjjriog^;!; 




[{x-hsy + y^+z^]^ 


and determining /(.v) from the integral etpiation 

y,(ir) 


1 r' f{s)di< 


^^0 

which is obtained by putting y = z^O. The integral equation is solved most cojiveniently by 
using the well-kiiown expansion 


Pn(s)<2n(r) 


and tlie integral formula 


It is thus evident that 


Jj 


Pm(s)Pt.(s)ds = 


0 m±n 
2 m=^n 
2 n4 1 


J{s)^F,{s)=^s. 

The strength of the source associated with an element hds is 


Ah, sd^. 47Tp 


Multiplying this by Ux = Uhs and integrating with regard to .s between - 1 and 1 , we get 

4tt ,,, rj 4ir<) ah^U^ 

3 3 ^ 1 , 1+6 
1 - e* 2 e 1 - c 

The kinetic energy of the fluid plus the kinetic energy of the fluid displaced is, on the other 


hand 



3 2 L 2-OoJ 

and 

,, 2(l-e’)r 1 1, 1+el 


Thus Munk’s theorem is again confirmed. 
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In the case of a prolate spheroid moving broadside on we have 


^ 1 


and the relation between K and Ic is 


where I is the perimeter of the meridian section. The potential (^b lURy be expressed in the 
forms 


where 


where 


'2-% I (a2 + X)K6^-HX)^ 


a» + X“^6^ + X 




1 1 + e e(l-2e^) 

alogr-e- 1-e’ 


= -AV 


<^>b=- 




{l — s^)ds 


At Doctor Xluiik's suggestion one may interpret these results with the aid of the idea of 
complete momentum, i. e., the momentum of the fluid plus tha momentum which the fluid 
displaced by the body would have if it moved like a rigid body with the same velocity as the 
body. 

Let .¥a and ¥b denote the complete masses for motions parallel to the axes of x and y 
respectively, then 

4 4 

¥a = girpa&c (1+W Mh= ^TTpahc (1-f^b) 


and we may write 


3¥a^ sds 

Sirh [ (x - As) ^ H- 


SMbVy 


n+i 

' I ^ Jl-®' 
J ^[(x-As)* + 


— s^)_ds 


These equations show that when the Complete inomentum is given the velocity potential $ 
and the sources from which it arises are the same for a series of confocal spheroids.® This is 
true for any angle of attack as is seen by superposition. This revolt is easily extended to the 
ellipsoid, for we may write 


where 




d\ 

(a* + X) a' 


i ^ _ 

a> + Xo'^6^+Xo 


« This is an extension to throe dimensioni of a theorem that has been proved for the elliptic cylinder. Cf. Max. M. Munk, Notes on Aero- 
dynamic Forces. T^hnical Note No. lOt, National Advisory Committee for Aeronautics. 
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It is easily seen that T is the same for a system of confocal ellipsoids. This result may be 
used to find an appropriate system of singularities distributed over the region bounded by a 
real confocal ellipse, the result is the same as that already found for the elliptic disk.^ 

It is well known that an ellipsoid has three focal conics, one of which is imaginary, and the 
question arises whether there is more than one simple distribution of singularities which will 
produce the potential. This question will be discussed in Section III. 

When a prolate spheroid is spinning with angular velocity Qy about the axis of y, the 
velocity potential 4> is given by the formulae 

= (l-M’)i log 

*/-l 

where .A is a constant to be determined by means of the boundary condition 


Sin w= ~ r. Az 




U-s’) ds 




It is easily seen that 


54 > 


^ dz\ 


A ^ (2 — e’) log i , * ,1 = aVil, 

^ ^ ^l-e e l-e^j ^ 


The energy may be expressed in terms of the mass of the fluid displaced by means of the 
formula 


(the coefficient having been tabulated by Lamb) or it may be expressed in terms of other 
characteristics with the aid of our coefficient The values of the various coefficients k and 
K are given in Table I. The suffixes a and c are used to indicate the axis along which the 
spheriod is moving. The coefficients h and refer to the case of rotation. It will be seen 
that the coefficients K vary only slowly and the same remark applies to the product (/ +ia) 
(/ _}_ ^ One advantage in using the coefficients Ks. and Kc is that it is not necessary to compute 

the volume of the hull of the airship. Since Kc varies very slowly indeed when the fineness 
ratio a/c is in the neighborhood of 6, it follows that if we take Kc = 0,6 for an airship hull we shall 
not be far wrong. 

Table I. 


a 

c 


1 

ke 

(1+*.) 

(l+*o) 

0.500 

Xc 

0.500 

k ' 


1.00 
1 155 

0.500 

0.500 

2,250 

0 

1 

0.906 1 

1.50 

0.305 

0.821 

2.116 

0. 457 

0.523 

0.244 


2.00 

0.209 

0.702 

2.058 

0.418 

0.541 

0.399 

0.895 

2.005 


! 





0.685 

2.09 

0.122 

6.^ I 

2.023 

0.3 W 

0.571 

0.582 


3. 571 




...1 1 

1 


0.512 

3.99 

0.082 

0.880 

2. 012 

0.317 1 

0-587 1 

0.689 


4.99 

0.059 

0.895 

2.006 

0.294 

0.599 



, 6.01 

0.045 

0. 918 

2.001 

0.270 

aeo 6 

0.807 


6. 97 

0.038 



0.250 

1 * 



8. 01 

0.029 



0.232 

i 

9.02 

0.024 



0. 216 


9. 97 

0. 021 


0.^ 


1 

1 

1 ^ 

1 

2 

0 

0. 637 

1 

0. 477 


In this table use has been made of the coefficients computed by Lamb, It should be 
noticed that £. + 2^0 is very nearly constant for values of ^ lying between 1 and 6. This fact 
may be used to compute Kc.y^hen is known using a formula such as 

£. = •743 £. 

The value thus found is too large for large values of ^ and too small for small values of ^ • 


1 Cf. Lamb's HydrodjTiamics, 3d ed., ch V, p. 145. 
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3. Oblate spheroid . — In the case of an oblate spheroid moving with velocity U in the 
direction of its axis of symmetry, which we take as axis of x, we have 


2r, • , 1 


where e is the eccentricity of the meridian section. In this case 

b = c, a = c^JT^ e^. 

The velocity potential is 


where 


^<x> 

_ac^Ux I 

Ju 


d\ 

^^ + X)*(c^-hX) 


a^-f-Xp'^c^ + Xo 


Introducing the spheroidal coordinates 

x = h^i^y y = (ji cos Wy z = (ji sin w, oy = h + 1)^ h — ce. 

(1-f cot-if) 


we may write 
where 


a^lc^ — a? 


~cos’^ - \= —h^ U. 


We also have 


f r % 

Jq 4/0 


= {y — hs cos wY-^ {z~-hs sin wY-\-x^. 

When an oblate spheroid moves with velocity TF at right angles to its axis of symmetry 
we have 


/3o = 7o = 


3 ^ [sin"* « — «Vl 


and the relation between tc and is now 

TT- IdCfy 

^^~2Ta 

The velocity potential <f>o is given by the formulae 


4>c = \ 


f 


d\ 


2-7, I (a’ + X)4(c^ + XF 


= ^(1 -M^)* (f*+ 1)* M 


= a[2sVI - s^-dsf |(i) 

Jo Jo 

R^={y-h8 COS {z — hs sin + h^ce^ — 

A |cos-“^ ^ ^ Vc* ” ^1*1 — 
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Some values of the coefficients fcc, Koy are given in Table II. 

TABLE II. 


e 

a 


AT. 

1.00 

0.500 

0.500 

1.50 

0.621 

0.523 

2.00 

0.702 

0.541 

2.99 

0.803 i 

0.571 

3.99 

0.860 

0.587 

4.99 

0. 895 1 

0.599 

6.01 

0. 918 

0.606 

00 

1,000 : 

0.637 



Ko 


0.500 

0.500 

2.250 

0. 382 

0. 4824 

2.240 

0.310 

0. 477 

2.229 

0.223 

0. 473 

2.205 

0. 174 

0. 474 

Z184 

0. 143 

0. 477 

2. 166 

0, 1205 i 

0. 478 

, 2. 149 

0.000 i 

0.500 

2.000 


When an oblate spheroid spins with angular velocity about the axis of y the velocity 
potential <f>^ is given by the formulae 


We also have 


(j>^= Czx 


f 


d\ 

(c^ + X) (a^+X)A 


(c^ — a^)^abco>y 

“ 2 (c^ - a^) + (c^ 4- a^) (yq ^ a©) 

(1 - m’) * (f" + 1) ‘ [3r cot->f - 3 + J 


A\S^ .,€ 08 -^-- 


c—a^ 


c c' -yj^ - a? 




4>'=- 


«/0 4/0 


III. THE METHOD BASED ON THE USE OF SOURCES AND SINKS. 

It was shown by Stokes * that the velocity potential for the irrotational motion of an 
incompressible nonviscous fluid in the space outside a sphere of radius^ a, moving with velocity 
Uj is the same as that of a doublet of moment 2-k Ua^ situated at the center of the sphere. This 
result has been generalized by Rankine,® D. W. Taylor/^ Fuhrmann/^ Munk/^ and others, two 
sources of opposite signs at a finite distance apart giving stream lines shaped like an airship. 

Munk has shown in a recent report that the intensity of the point source near one end of 
an airship hull may be taken to be r^ir U, where r is the radius of the greatest section of the 
ship and ir the distance of the point source from the liead of the ship. The total energy of 
the fluid displaced is then 

r=g7rr’pZ7^ 

and the apparent increment of mass of the airship is equivalent to about 2^ per cent of the 
mass of fluid displaced. 

In this investigation the airship is treated as symmetrical fore and aft, the two sources 
of opposite signs being equidistant from the two ends and the contributions of the two sources 
to the kinetic energy being equal. The, final result is identical with that for an elongated 
spheroid with a ratio of axes equal to 9. 

• Cambr. PhU Trans., vol. 8 (1843). {Math, and Phys. Papers, Vol, I. p. 17.{ 

> PhU. Trans. London (1871), p. 267. 

w Trans. British last. Naval Architects, voi. 35 (1894), p. 385. 

Ji Jahrb. der Motorluftschlfl-Studiengesellschaft, 1911-12. 

National Advisory Committee for Aeronautics, Reports 114 and 117 (1921). 
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It is thought that a lack of fore and aft symmetry will still further reduce the values of the 
coefficients Ic and K, To get an idea of the effect of a lack of symmetry we shall consider the 
case of a solid bounded by portions of two orthogonal spheres. In this case, as is well known, 
the velocity potential may be derived from three collinear sources. We may in fact write 


^ = 2 a’ ^ 


. JJ 

r 2 


V 

r 2 


cos 

sin^ 


2?’ V 

-\tu 


COS 0 

W 

sin^ 0 
R 


where a and a' are the radii of the two spheres (r, 6), (r' {R, 0) are polar coordinates re- 

ferred to the three sources as poles, the angles being measured from the line joining the three 
sources. If $ is a common point of the two spheres R is measured from the foot of the per- 
pendicular from Q on the line of centers, while r and r' are measured from the centers of the two 
spheres respectively. The quantity p represents the distance of Q from the line of centers and is 
given by the equation 


By means of Munk's theorem we infer that the complete energy is given by the formula 
2T=2irp {l^]c)VU^ 

where 


V=l^2(a^+ +2a» +2a'» - 3 

is the volume of the fluid displaced. 

The fineness ratio, i. e., the ratio of the length to the greatest breadth, is 

y. d ~j-(t 

* 2a 

Some values of h and K are given in Table III and curves have been drawn in Figure 2 to 
show the effect of a lack of fore and aft symmetry. For a comparison we have given in Table 
III the values of h and K for a spheroid of the same fineness ratio. The high value of K for 
the two orthogonal spheres is undoubtedly due to the presence of a narrow waist. The sudden 
drop in the value of K is probably due to the lack of fore and aft symmetry. The coefficient 
K shows the effect of a change in shape much more clearly than the coefficient Ic, 

TABLE III. 


a 

k 

K 

/ 

k (spheroid). 

iCCspheroid). 

1.0 

0.313 

0. 5897 

1.707 

0. 243 

0. 440 

0,9 

0.315 

0.5136 

1.622 



0.8 

0.329 

0. 4708 

1.54 



0.75 

0. 334 

0.4509 i 

1.5 

0.305 

0. 457 

0.66 

0.363 

0.4603 

1.434 



0.41 

0. 448 

0. 471 

1.25 



0.29 

0.48 

0.488 1 

1. 166 



0 

0.5 

0.5 ! 

1 

1 

0.5 

0.5 


It appears from an examination of the case of the oblate spheroid that the motion of air 
round a moving surface of revolution can not always be derived from a number of sources at 
real points on the axis. For the oblate spheroid the sources, or rather doublets, are in the 
equitorial plane. It is possible, however, to replace these doublets by doublets at imaginary 
points on the axis as the following analysis will show. 

23—24 n 
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if F{x, z) is a potential function, we have the equation 


n2it 


y — o cos o),x — (T sin wj aw 


■■p 


[a; + i(T cos X,y,z] dX 


which holds under fairly general conditions. On account of this equation we may write 

p2r n2ir 

I/O t/O c/0 t/O 

s f /.")* f id) f f !(«-)"■ 

t/o c/0 c/0 c/0 


where 
Putting 
making the substitution 


.1 f g d) d<r r Q) = gd)d A Q,) 

c/0 c/0 c/o t/0 




(y — o' cos o))^-f (2 — <7 sin w)^, (/-f- icr cos + + 

/ (^) =/i (<^) = ^ ^ •7'') h y (c^) = O' ^ 

(7COsA" = {, dXVo"^ — 
changing the order of integration and making use of the equations 


i ‘ 

r 




17,1 


16^ 

which are easily verified by means of the substitution 

cos" ^ + sin" e, 

we find that the potentials for the oblate spheroid in the three types of motion may be written 
in the forms 

A 


= “ ~ G") 

J' di = + 2h^J^J “■ dz (g") 


“ 8 F 


where 


fl"-[(x + i|)' + y^ + 2*]h 


These formulae resemble those for the prolate spheroid. 

A distribution of sources or doublets over the elliptic area bounded by a focal ellipse of an 
ellipsoid may be replaced by a system of sources or doublets at imaginary points in one of the 
other planes of symmetry by making use of the equation.'* 

*• H. Bateman, Amer. Journ. of Mathematics, vol. 34 (1912), p. 335. 
w H. Bateman, loc. clt., p. 336. 
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F[x~<j COS B COS a, sin 6^ z] dS^ 



F[x+ i(T sin X sin a, 0 + ?V cos X] dX 


which likewise holds under fairly general conditions when F{Xf y, z) is a potential function 
and a an arbitrary constant. 

The theorem relating to the transformation of doublets in a central plane into a series of 
doublets at imaginary points on the axis of symmetry may be written in the general form 



d d d 
bi by bz 



^ b d\ 1 
dx by bz) ir^ 


F (i) dt 


wliere the functions/(cr) and F{^) are connected by the integral equation 


F{k) = 



ficr) da 


In order that Munk’s theorem may be applicable to doublets at imaginary points as wel 
as to doublets at real points we must have the equation 


Now 



h 

f {a) da 


hence the formula is verified and the complete mass may be calculated from doublets at imaginary 
points by adding the moments and using Munk^s formula. 


IV. CASES IN WHICH THE MASS CAN BE FOUND WITH THE AID OF SPECIAL HARMONIC FUNCTIONS 

It is known that the potential problem may be solved in certain cases by using series of 
spheroidal toroidal, bipolar, or cylindrical harmonics. Thus it may be solved for the spherical 
bowl, anchor ring, two spheres,*^ and for the body formed by the revolution of a lima(;‘on about 
its axis of symmetry. The last case is of some interest, as it indicates the eflect of a flattening 
of the nose of an airship hull. Writing the equation of the limacon in the form 


r = 2a? 


s + cos 6 
~ s^~l 


where r and B are polar coordinates, we find on making the substitutions 

r cosB = x=^^ — r)^. r sin S = y = 2^T] 

a sinh a ^ X 

^ cosh a — cos X ^ cosh a — cos x 

that the potential for motion parallel to the axis of symmetry is 


<t> = 




[F m (cosh a) ^ m+1 (cos x) ^m+1 (cosh <r) Pm (cos x)] 

I’ For referoQces scj» Lamb’s Hydrodynamics, 3d ed., pp. 126, 119; and A. B. Basset, Hydrodynamics, Cambridge, 18SS, Vol. I. 
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where Pm W and Qm (s) are the two types of Legendre functions (zonal harmonics) and P'm 
(s) are the derivatives of Pm (s) and Qm is) respectively. The stream-line function as found 
by Basset is in our notation. 




• V, ( 2 m + 3 ) 

cosh <r — cos X 


^'m+i (s) 


P'm^i (cosh o) P'm^l (cos x) sinh^ a sin’ x 


At a great distance from the origin we have the approximate expressions 




P ’ 




U 






R 


= cosh a — cos x> 


oo 



( 2 m + 3 ) (m+ 1 )^ 


<2 m., 

P'm + 1 


( 5 ) 

(s) 


which give the sum of the moments of the doublets from which the potential arises. The 
coefficients A: amd K may now be calculated with the aid of Doctor Munk's theorem and an 
incomplete table of spheroidal harmonics which is in the autlior’s possession. We thus obtain 
the values 

TABLE IV. 


s 

’ / 

k 

K 

fc (spheroid) K (spheroid) 

00 

I 1 

1 

0.500 

0.500 

0.500 

0.500 

3 

1.05 

0. 527 

0. 507 

0.5i2 

0.502 

2 

1.10 

0. 54H 

0.513 

0. 521 

0.505 

1.2 

i 1. 153 

0. 569 

0.518 

0. 530 

0.507 

1.1 

1. 154 

1 0. 573 

0. 523 

0. 516 

0. 507 

1 

1. 155 

0, 578 

0. 527 

0. 536 

0.507 



The corresponding values for an oblate spheroid are given for com- 
parison. The case in which s = .^ is particularly interesting because 
the limacon then has a point of undulation at the nose. When s<2 
the limacon curves inward at the front, as may be seen from the dia- 
grams in Figure 3, and the apparent mass is probably increased on 
account of fluid being confined in the hollow. In calculating the fine- 
ness ratio in such a case the length has been measured from the rear 
to the point where the double tangent meets the axis. 


The values for the spheroid have been obtained by interpolation from Table II. The values of k and K for thec ardioid s- 1 have been esti- 
mated by extrapolation. 
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